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ABSTRACT 

We have discovered the pulsar associated with the supernova remnant G21.5-0.9. PSR J1833— 1034, with 
spin period P = 61.8 ms and dispersion measure 169cm~ 3 pc, is very faint, with pulse-averaged flux density 
of w 70 /iJy at a frequency of 1 .4 GHz, and was first detected in a deep search with the Parkes telescope. 
Subsequent observations with Parkes and the Green Bank Telescope have confirmed this detection, and yield 
a period derivative P = 2.02 x 10~ 13 . These spin parameters imply a characteristic age t c = 4.8 kyr and a spin- 
down luminosity E = 3.3 x 10 37 ergs s" 1 , the latter value exceeded only by the Crab pulsar among the rotation- 
powered pulsars known in our Galaxy. The pulsar has an unusually steep radio spectrum in the 0.8-2.0 GHz 
range, with power law index ss 3.0, and a narrow single -peaked pulse profile with full- width at half maximum 
of 0.04P. We have analyzed 350 ks of archival Chandra X-ray Observatory High Resolution Camera (HRC) 
data, and find a point-like source of luminosity « 3 x 1Q~ 5 E, offset from the center of an elliptical region of 
size w 7" x 5" and luminosity w 10~ 3 .E within which likely lies the pulsar wind termination shock. We have 
searched for X-ray pulsations in a 30 ks HRC observation without success, deriving a pulsed fraction upper 
limit for a sinusoidal pulse shape of about 70% of the pulsar flux. We revisit the distance to G2 1.5-0. 9 based 
on H I and CO observations, arguing that it is 4.7 ± 0.4kpc. We use existing X-ray and radio observations of 
the pulsar wind nebula, along with the measured properties of its engine and a recent detection of the supernova 
remnant shell, to argue that G21.5-0.9 and PSR J1833— 1034 are much younger than r c , and likely their true 
age is < 1000 yr. In that case, the initial spin period of the pulsar was > 55 ms. 
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1. INTRODUCTION 

"Crab-like" remnants of supernovae (SNe) display no read- 
ily observable manifestation of an ejecta-driven blast wave; 
instead, the relativistic wind of particles and magnetic fields 
released by a young pulsar is confined by the ambient pres- 
sure giving rise to a non-thermal compact nebula. Such pul- 
sar wind nebulae (PWNe), particularly with knowledge of the 
spin parameters of their central engines, are magnificent labo- 
ratories for the study of astrophysical shocks, the properties of 
pulsar winds, and the spin evolution and birth parameters of 
neutron stars. Eventually, depending on explosion energy and 
environment, high velocity SN ejecta are expected to generate 
a brightening shell surrounding the PWN, creating a "com- 
posite" supernova remnant (SNR). 

G2 1. 5-0.9, known for 35 years llAltenhoff et alJ 
119701 IWilson & Altenhofl 1 19701) and long-c lassified 
ftVilson & WeileJ 119761 iBecker & Szvmkowiakl 119811) 
as one of about ten Crab-like remnants dGreenl 120041) . is 
a bright and absor bed (neutral hydrogen column density 
N H = 2 x 10 22 cm" 2 ; ISafi-Harb et afll2001l) centrally peaked 
X-ray and polarized radio source. H I absorption measure- 
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ments (Caswell et al. 1975; Davelaar, Smith, & Becker 1986) 
place the system at > 4.4 kpc. Its radio and X-ray lumi- 
nosities are ~ 10 and ~ 100 times below thos e of the Crab 
Nebula, respectively (Helfan d~& Beckerlfl987l) . A break in 
the observed spectrum of G2 1.5 -0.9 was ori ginally reported 
at a frequency v b — 50 GHz (e.g.. lSalter et alJl989l) . although 

GHz 
120011 

X-ray observations in the modern era have advanced our 
knowledge of G2 1.5-0. 9 enormously. At the c enter of the 
system is a bright compact source (« 2" in radius; Sla ne et all 
2000) that may include the pulsar and delineate its wind ter- 
mination shock. Surrounding this is the previously known 
PWN (« 40" in radius), which is more compact than the 
radio nebula, and wi th a spectrum that steepens ra dially 
from the center ( e.g.. ISafi-Harb et alJl200lt IWarwick et all 
120011 iBock et alJ 120011) . But most remarkably, surround- 
ing the PWN is a low surface-brightness non-th ermal halo, 
» 15 0" in radius, with no radio counterpart (Sla ne et all 
2000), and apparently with a spa tially invariant photon in- 
dex ([Matheson & Safi-Harb 2005). The v ery dee p Chandra 
image presented by Matheson & Safi-Harb (2005) beautifully 
demonstrates evidence for limb brightening in this halo. It 
now appears likely that a portion of the emission observed 
in the halo originates in the long sought SNR shell, while it 
is possible that dust scatte ring also contributes significantly 
to the observed emission (Bocchi no et alJ l2005). especially 
within about 80" of the center. While unusual, such a non- 
thermal shell is n ot unprecedented (e.g., Kovam aet alJI 19951 
iSlane et alii 19991) . and could be due to relativistic electrons 
accelerated at the SN ejecta-ambient gas shock front. G21.5- 
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0.9 has also been detected in hard X-rays by INTEGRAL, with 
a 20-100 keV flux of ~ 3 millicra bs and a flatter spectrum 
than that of the Crab Nebula l lBird et all2004l) . 

Meanwhile, the pulsar powering the PWN at the cen- 
ter of the composite SNR G2 1.5-0. 9 ha s remained unde- 
tected, despite searc hes in the radio (e.g., Kaspi et alj [12261 
I Crawford etall 120021) and X-ray bands JSlane et alJ 12000: 
ISafi-Harb et all2001tljTP alombar alTMereghettil2002l) . We 
have now, i n the course of ext ensive deep radio searches of 
PWNe (e.g.. lCamilo et alJl2002alcT) . discovered PSR J1833- 
1034, the central engine of G21.5-0.9, which we report on 
here. After submitting this paper, we became aware of a n 
independent discovery of this pulsar by Gupta et al. (2005), 
whose results on the pulsed signal from this source largely 
agree with ours. 

We begi n by describing the pulsar search at radio wave- 
lengths (§ 12. 1> and its confirmation and follow-up observa- 
tions f§ 12.211 . We then present an X-ray imagi ng a nalysi s of 
G2 1.5-0.9 (§E3 and search for pulsations (§|fij. In §0 
we obtain the distance to the system, after which we discuss 
th e wi nd terminat ion shock region, nebular magnetic field 
(§ 13. 2> . and age (§ 13. 3i of G2 1.5-0.9. We compare the prop- 
erties of PSR J1833-1034 and its SNR to those of SNR 3C 58 
and its pulsar in § 13.41 and in § 13.51 we remark on the ra- 
dio spectrum of PSR J1833-1034, putting it in the context of 
those of other young pulsars. We summarize our work in §|4] 

2. OBSERVATIONS AND ANALYSIS 
2.1. Radio Pulsar Search 

We observed the position R.A. (J2000) 18 h 33 m 33. s 8, Decl. 
(J2000) -10° 34' 10", centered on G2 1.5-0.9, at the Parkes 
64-m radio telescope on 2002 January 6 for a full transit of 
8.0 hr with the central beam of the multibeam receiver system 
JStavelev-Smith et alll996l) . The observation was centered at 
a frequency of 1374MHz, with a bandwidth of 288 MHz for 
each of two polarizations split into 96 contiguous frequency 
channels. Total-power signals were detected and sampled 
with 1-bit precision every 0.8 ms, with a total of 36 million 
time samples for each of 96 frequency channels recorded to 
magnetic tape for analysis. 

We analyzed the d ata using standard pulsa r search 
techniques (see, e.g., lLorimer & Kramer! |20 i i 5j) imple- 
mented in the PRESTO software package (Ransom 1200 It 
[ Ransom. Eikenberrv. & M iddleditch 2002a). While the 
Cordes & Lazio ( 120021) model for the Galactic distribution 
of free electrons predicts a dispersion measure (DM) > 
282cm" 3 pc for the distance to G2 1. 5-0.9 of > 4.4 kpc, 
such predictions can be inaccurate (e.g.. lHalpern et alll2001t 
Camil o"et alJ |2Q02d), and we searched the DM range 0- 
2674 cm" 3 pc, up to twice the maximum expected Galactic 
value in this direction. We first analyzed the non-dispersed 
(i.e., DM = cm" 3 pc) data to identify strong interference sig- 
nals that we subsequently masked (in the time, radio fre- 
quency, or fluctuation frequency domains). We de-dispersed, 
barycentered, and searched the data for each of 524 trial DMs 
while maintaining close to optimal time resolution through- 
out the DM range searched. A pulsar with very large spin 
period derivative, P, may change its spin frequency percepti- 
bly within such a long observation, diminishing the sensitivity 
unless a number of trial acceleration values are used to gen- 
erate multiple time series to be searched: e.g., the Crab pul- 
sar's fundamental would move by 0.3 Fourier bins within our 
observation, and a young pulsar in a binary system could be 



TABLE 1 

Radio Observations of PSR J1833-1034. 



Date 


Telescope 


Time 


Frequency 


Bandwidth 


Backend" 


J 1 b 

J sys 


SpSR 1 " 


(MJD) 




(hr) 


(MHz) 


(MHz) 




(K) 


QjJy) 


52280 


Parkes 


8.0 


1374 


288 


filter bank 


40 


66 


53474 


GBT 


0.5 


1950 


600 


SPIGOT 


39 




53479 


Parkes 


5.9 


1374 


288 


filter bank 


40 


76 


53482 


GBT 


3.3 


820 


48 


BCPM 


83 


277 


53486 


GBT 


4.5 


1950 


600 


SPIGOT 


39 


13 


53502 


GBT 


0.4 


820 


48 


BCPM 


83 


283 


53525 


GBT 


1.0 


820 


48 


BCPM 


83 


181 


53557 


GBT 


0.8 


820 


48 


BCPM 


83 


193 


53575 


GBT 


2.8 


820 


48 


BCPM 


83 


171 



"For details on observing equipment, see: filter bank i Manchester et al. 2001 ); SPIGOT 
IRansom et all l2Tjo"5l Kaplan et al] f%)U3l : BCPM IBacker et all 119971 ICamilo et alJ 

b Inc ludes contribution from Galaxy (scaled from the 408 MHz map of Ha-laiiujLjilJ 
1982, with spectral index of 2.6) and SNR (« 6 Jy at 1 GHz with flat spectrum lS^fterTtTl] 

Tntegrated period-averaged pulse profile, compared to the rms of the off-pulse region, 
scaled via the radiometer equation using the observation bandwidth and time, and assum- 
ing the nominal telescope gain and system temperature (previous column). For SPIGOT, 
we used only 1650-2050 MHz data to calculate flux density (see caption to Fig.QJ. These 
estimates have an uncertainty of about 25%, exclusive of variations due to interstellar 
scintillation. 



affected even more by orbital Doppler accelerations. For this 
reason we used Fourier-domain acceleration searches to cover 
a range of P that would cause a signal to drift by up to ±16 
Fourier bins for the highest of 1, 2, 4, 8, and 16 harmonics. 
We summed these harmonics incoherently to retain sensitivity 
to narrow pulse profiles. 

We detected a clear pulsar signal at DM rs 170 cm" 3 pc with 
barycentric period P = 61.84 ms, P < 5 x 10~ 13 and a narrow 
profile with full-width at half maximum of 0.05P (Fig. [2 mid- 
dle panel). 

2.2. Confirmation and Follow-up Observations 

On 2005 April 14 we observed the pulsar with the 100-m 
Green Bank Telescope (GBT) for 0.5 hr at a central frequency 
of 1950 MHz. While the pulsar was not detected in an unbi- 
ased way in this short observation, searching a range of pe- 
riods for assumed values of P < 5 x 10~ 13 yielded the most 
promising candidate for P ~ 2 x 10" 13 . On 2005 April 19, 
we clearly re-detected the pulsar at 1374 MHz at Parkes with 
P = 61.86ms, implying P = 2.02 x 10~ 13 . The retrospective 
detection at GBT also implies that its positional uncertainty 
from radio observations is about ±3' (telescope beamwidth). 

We have now observed the new pulsar at 820 MHz, in ad- 
dition to the two frequencies described above (see Table Q. 
There is some evidence for fluctuation in received flux den- 
sity, likely due to interstellar scintillation (see Fig. [Q. The 
flux densities listed in Table ^ are calibrated only approxi- 
mately, and are therefore subject to some future adjustment. 
Nevertheless, we infer that the spectrum of the pulsar is steep, 
with a weighted fit to the values presented in Tabled giving 
a = 3.0 ±0.3 (throughout this paper, S„ oc v~ a ). The flux den- 
sity of ~ 65 mJy at a frequency of 610 MHz estimated by 
Gupta et al. (2005) is consistent with this determination of a. 

The new pulsar has a magnetic field strength at its 
light cylinder of B\ c = 1.4 x 10 5 G, the sixth highest value 
among pulsars kno wn in the Galactic disk, ex aequo with 
PSR J2229+6114 (lHalpern et alJl200l . Since the emission 
of "giant" radio pulses appears to be correlated with B\ c (e.g., 



PSR J1833-1034in SNR G21.5-0.9 
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TABLE 2 

Measured and Derived Parameters for PSR J 1 833—1034. 



Parameter 


Value 


Right ascension, R.A. (J2000) 


18 h 33 ra 33?57 


Declination, Decl. (J2000) 


-10°34'07;'5 


Spin period, P (s) 


0.061865672117(5) 


Epoch (MJD) 


53545.0 


Period derivative, P a 


2.02025(3) x lfr 13 


Period derivative, P b 


(2.0206 ± 0.0004) x 1CT 13 


Dispersion measure, DM (em -3 pc) 


169.5(1) 


Post-fit rms timing residuals (ms) 


0.23 


Data span (MJD) 


53479-53575 


Pulse full-width at half-maximum (ms) 


Ri2.5 


Flux density at 820 MHz (jjjy) 


^221 


Flux density at 1374 MHz, Shoo (/Ldy) 


»71 


Flux density at 1850 MHz QiJy) 


13 


Spin-down luminosity, E (ergs s _I ) 


3.37 x 10 37 


Surface dipole magnetic field strength (gauss) c . 


3.58 x 10 12 


Characteristic age, t c (kyr) 


4.85 


Distance, d (kpc) 


4.7 ±0.4 


Radio luminosity at 1.4GHz, Snoad 1 (mJykpc 2 ) 


« 1.5 


Radio spectral index, a 


« 3.0 



NOTE. — Celestial coordinates are assumed to corresp ond t o the central peak 
of the X-ray PWN seen in Chandra HRC observations (§ 12.31 Fig.|5}. Numbers 
in parentheses are three times the nominal 1 a TEMPO uncertainties in the least- 
significant digits quoted. 



"This P was determined from the timing solution based on data collected during 
the span listed. 

b This P was obtained from the TEMPO-derived P for MJD 53545.0 and the 
period measured from the discovery observation on MJD 52280. 
C B P = 3.2 X 10 19 (P^) 1/2 G (see, e.g. JManchester & TavloJl577l> . 

ICognard et al.lH996t iJohnston & RomanH l200l . we made a 
preliminary search for such pulses in the long 820 MHz data 
set (Tabled, but have not detected any. 

We measured pulse times-of-arrival (TOAs) from the GBT 
and Parkes data and fitted a simple spin-down model to the 
2005 TOAs using TEMPO 8 . The results of the fit are shown in 
Table|2] The position used for the pulsar was that determined 
from Chandra observations (§ 12.31 . The P measured from 
the phase-coherent solution at the present epoch matches that 
derived from two period measurements separated by 3.5 yr. 

2.3. X-ray Imaging of G21. 5-0.9 

G2 1.5-0. 9 is used as a calibration target for the Chandra 
X-ray Observatory, and multiple observations exist using a 
variety of detector modes and aimpoint positions. For a set of 
six observations carried out with the Advanced CCD Imag- 
ing Spectrometer (ACIS), with the target placed within 1' of 
the optical axis of the telescope and falling on CCD S3, we 
used the CIAO 9 task merge_all to generate a combined im- 
age of the PWN (Fig. |2j. The integrated exposure for the 
image is 58 ks, and we have applied a Gaussian smoothing 
of two pixels (0"98). The image shows an extended core, 
identi fied in earlier C handra observations ( Slane et al J l200Ct 
ISafi-Harb et all 2001). surrounded by a broader nebula with a 
radius of « 40". There is considerable structure in the outer 
nebula region, including filamentary features and a distinct 
low-surface-brightn ess "notch" in the outer northwest of the 
PWN (see alsolSafi-Harb et alJl200ltlMatheson & Safi-Harbl 
120051) . We split the image in Figure [2] into two energy bands 
(0.3-6.0keV and 6.0-10.0keV) and computed radial surface 

8 See http://pulsar.princeton.edu/tempo 

9 See http://cxc.harvard.edu/ciao 



brightness profiles: the profile in the hard band is slightly 
steeper than that in the soft band, although X-ray emission 
extends in both bands to roughly the same outer radius. 

To further investigate the compact core, we created a 
merged image from 21 observations with the Chandra High 
Resolution Camera (HRC), using data from both the imaging 
(HRC-I) and spectroscopy (HRC-S) cameras (the latter used 
here without transmission gratings). The emission-line star 
SS 397 was identified in the field of each observation, and its 
2MASS position (R.A. [J2000] 18 h 33 m 27?70, Decl. [J2000] 
-10°35'23"0) was used to register separately each individual 
observation in order to improve the overall astrometry. The 
resulting merged image of the core region, corresponding to a 
total integration of 35 1 ks, is shown in Figure[3] With a typical 
2MASS positional uncertainty of < 0"2 and determination of 
the X-ray centroid to better than an HRC pixel, we estimate 
that any residual blurring/registration error in this image is 
at the l evel of < 0"3. A com parison of Figure [3] with Fig- 
ure 4 of lSafi-Harb eTail ll200"l reveals significant differences 
in morphol ogy, plus an off s et of ~ 1 "5 in the position of peak 
emission. Safi-Harb et al J J2001I) did not discuss any efforts 
to correctly register each exposure, and it is possible that their 
image has been affected by slight shifts between observations. 

Immediately evident in Figure [3] is a compact source sur- 
rounded by an elliptical emission region of size w 7" x 5" 
(in a very saturated image, this core appears to be about 50% 
larger in each dimension; see, e.g., Fig.|2li. This is to be com- 
pared with the 90% encircled-energy radius of ss 1" for HRC. 
The point-like source is located at R.A. (J2000) 18 h 33 m 33?57, 
Decl. (J2000) -10°34'07."5 and offset from the center of the 
ellipse, possibly indicating a Doppler-boosted toroidal struc- 
ture surrounding a pulsar, as is observ ed for the Crab and 
other pulsars imaged with Chandra (e .g.jN g & Romanil2004t 
iHessels etail2 004: Slan e et alJ200 4). In this scenario, the in- 
ferred pulsar spin axis projected on the sky is perpendicular to 
the major axis of the ellipse, at a position angle of ~ 40 north 
throug h east (sim ilar to the position angle first suggested by 
Fiirst et al. 1988, for the central source in G21. 5-0.9). The 
major axis of the ellipse is aligned toward the notch of re- 
duced e mission on the northwest rim of the PWN (Fig. |2j 
see also Matheson & Safi-Harb 2005). We note the similar- 
ity with the "bays" observed in the Crab Nebula, which are 
rough ly aligned with the tor us surrounding its pulsar (e.g., 
iFesen. Martin. & Shulil9"92l) . 

In order to estimate the flux of what is presumed to be the 
pulsar, we used the CIAO package Sherpa to perform spatial 
modeling in which we fit the core region of G21 .5-0.9 (Fig.|3j 
to a model consisting of a point source (convolved with the 
PSF, calculated using the CIAO task mkpsf for an energy of 
1 keV), a 2-D Gaussian, and a constant background. While the 
extended emission is clearly not a 2-D Gaussian (we expect 
perhaps a Doppler-brightened torus superimposed on more 
uniform emission), this is the most suitable model available 
in Sherpa and provides a good fit. As expected, the best-fit 
values for the point-source position and Gaussian centroid are 
very different, limiting the point source flux assigned to the 
Gaussian component. In any case, there is a systematic uncer- 
tainty in this flux determination in excess of the formal fit un- 
certainty of 10%, and a more realistic geometrical modeling 
is required to improve significantly on this result. In the en- 
ergy range 0.5-10keV, the resulting count rate for the pulsar 
is 1.5 x 10~ 3 s , corresponding to an unabsorbed flux of about 
2.8 x 10~ 13 ergscm" 2 s" 1 for a spectral index of ay = 0.5, as- 
suming N H = 2.2 x 10 22 cm -2 ISafi-Harb et alJl2001l) . or a lu- 
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minos ity of about 9 x 10 32 ergs s -1 at a distance of 5 kpc (see 
§ 13. U . The assumed Nh may be somewhat in error, since sim- 
ple fits to the data yield values of Nh that differ by as much 
as 20% depending on radial position within the nebula. How- 
ever, such a variance in Nh introduces a variation of only 20% 
on the point source flux, smaller than the uncertainty arising 
from the difficulties involved in spatial modeling. Through 
such an absorbing column and at a distance of 5 kpc, thermal 
radiation from a lOkm-radius neutron star with a surface tem- 
perature of 1 .4 MK (a ppropriate for standard cooling to an age 
~ lOOOyr; see § I3.3> would contribute an HRC count rate of 
only ~ 10% that observed, and we do not consider it further 
here. 

The entire elliptical region of X-ray emission shown in 
Figure [5] including the compact component, has a 0.5- 
lOkeV count rate of 6.1 x 10~ 2 s _1 . Assuming ax = 1.0 and 
Nh = 2.2 x 10 22 cm" 2 , this corresponds to an unabsorbed 0.5- 
lOkeV flux of 1 x 10 -11 ergs cm" 2 s" 1 or a luminosity at 5 kpc 
of 3 x 10 34 ergss _I . 

2.4. Search for X-ray Pulsations 

G2 1.5-0. 9 was observed for 30 ks with the Chandra HRC- 
S on 2002 November 12 (ObsID 2755), with the detector 
operated in timing mode. In this mode, event arrival times 
are adjusted to compensate for a wiring error, allowing pho- 
tons to be time-tagged to a resolution of 16 /is. We cor- 
rected these event times to the Solar System barycenter us- 
in g the CIAO task axbary and the nominal pulsar position 
(§ 12.31 . We extracted 558 events from a circular region of 
radius 1", centered on the position of the compact source 
identified in the merged HRC image (Fig. |3j- We searched 
these events for pulsations around the spin frequency pre- 
dicted by the radio ephemeris (Table [2), covering approxi- 
mately 3 independent Fourier spacings (10~ 4 Hz). We used 
the Discrete Fourier Transform, which is optimally sensi- 
tive to sinusoidal pulsatio ns, and the Bayesian technique of 
iGregorv & Loredol ( 119921) . which is sensitive to pulsations 
with a more complex but unknown shape. In neither case 
did we find any evidence of X-ray pulsations. We estimate 
a 99% confidence-level upper limit on the pulsed fraction for 
sinusoidal pulsations of 21% , using the method described by 
iRansom. Gaensler. & Slanel l l2002bl) . If the pulsar has a sig- 
nificantly narrower X-ray pulse profile, as in the radio, the 
pulsed fraction would be significantly less than 21% (with the 
exact value depending sensitively on the pulse shape). This 
limit pertains to pulsations from the total observed flux within 
the extraction region and does not apply to the pulsar's flux 
alone. Since the pulsar likely contributes < 30% of the counts 
in the region considered, our upper limit represents a signifi- 
cant fraction (~ 70%) of the total signal from the pulsar, and 
hence is not especially constraining. 

3. DISCUSSION 

The P and P of the newly discovered pulsar imply a charac- 
teristic age t c = P/(2P) = 4.8 kyr and spin-down luminosity 
E = 4tt 2 IPP~ 3 = 3.3 x 10 37 ergss _1 , where the neutron star 
moment of inertia is / = 10 45 gcm 2 . T h is £ i s remarkably 
similar to that predicted by iSlane et aTl (1200 0) for the pul- 
sar in G2 1.5-0. 9 based on the X-ray energetics of the PWN, 
and is second only to that of the Crab among known Galac- 
tic rotation-powered pulsars. Still, until X-ray pulsations are 
detected from the core of the PWN with P = 61.8 ms or un- 
til the positional uncertainty of the pulsar is reduced to a few 



arcsec, it could be argued that it may not be the counterpart of 
the PWN. However, if the newly discovered pulsar were not 
associated with G2 1.5-0. 9, we would certainly expect to de- 
tect it in the deep Chandra observations of the field, given its 
relatively small inferred distance (§| 3.1> and lar ge E, and no 
such p otential counterpart exists (e.g., Matheson & Sa fi-Harbl 
2005). Therefore, this is clearly the neutron star responsible 
for the PWN in G21.5-0.9, which we designate PSR J1833- 
1034 based on the position of the central compact source de- 
tected with Chandra (Fig.[3J see also § !2.3l and Tabled- 

3.1. Distance to G2 1.5-0.9 and PSR J1833-1034 

To derive an upper limit on the distance to G2 1.5-0. 9, we 
consider the bright extragalactic source PMN J1832-1035, 
which lies just 18' from G2 1.5-0. 9. An H I absorption spec- 
trum of PMN J1832-1035, using data taken from th e Very 
Large Array Galactic Plane Survey dTavlor et alJ20 02). shows 
a strong absorption feature at a local standard of rest (LSR) 
veloci ty of +76kms~ 1 , but which is not seen toward G21.5- 
0.9 dCaswell et alJU975l Davelaar et al.lll986l) . The CO sur- 
vey of Dame, Hartmann, & Thaddeus (2001) shows a bright 
coherent region of molecular material at this velocity, present 
in the direction of both PMN J1832-1035 and G21.5-0.9. 
Thus, it is highly likely that there is cold neutral gas at an 
LSR velocity of +76kms _1 toward both sources. The lack of 
Hi absorption against G21.5-0.9 at this velocity implies that 
the SNR is in front of this cloud. 

We next compare th e emission spectrum of this CO to the 
H I emission survey of Hartm ann & Burtonl (119971) . In the di- 
rection of PMN J1832-1035, the H I emission shows a narrow 
feature of reduced emission in the spectrum at a velocity of 
+76kms _I . We interpret this feature as Hi self-absorption, 
indicating that cold gas on the near side of the tangent point is 
absorbing H I emission produced on the far side (Liszt, Bur- 
ton, & Bania 1981; Jackson et al. 2002) Usin g the Galac- 
tic rotation curve of Fich. Blitz. & Stark ( 1989) and assum- 
ing a Galactocentric radius for the Sun of 8.5 kpc, the dis- 
tance to this cloud and hence an upper limit on the distance 
to G2 1.5-0. 9 is 4.9 kpc. Combined with previous Hi mea- 
surements that put a l ower limit on the SN R's systemic LSR 
velocity of 65 km s -1 ( Davelaar et alll986l) . corresponding to 
a distance limit of 4.4 kpc, and adopting an uncerta inty in 
the systemic velocity of H i clouds of ±7 km s -1 dShaver et alJ 
H9821 iBelfort & Crovisieri rT984). we conclude that the kine- 
matic distance to SNR G2 1.5- 0.9 is 4.7 ± 0.4 kpc. 

The Cordes & Lazio (2002) model yields a distance d = 
3.3 kpc for PSR J 1833-1034, which is clearly an underestima- 
tion of the true distance, but supports a relatively small value. 
ISafi-Harb et alJ ( 120011) have estimated 4. 1< d < 5.2 kpc from 
a calibration to a range of X-ray-derived Nh values. Hereafter 
we scale the distance to PSR J1833-1034 and G21.5-0.9 in 
terms of d$ =d/{5 kpc). 

The X-ray spectra of the various components of G2 1.5-0. 9 
imply a foreground absorbin g column Nh ~ 2 x 10 22 cm" 2 
( War wick et alJ 1200 ll ISafi-Harb etafl 1200 ll) . Comparison 
with the DM given in Table implies a ratio Afa/DM » 40. 
This is substantially highe r than seen toward a ll but ~ 3 other 
pulsars (see discussion in Gaensle r et al.ll2004l) . Gaensler et 
al. argue that these high ratios indicate that a source is behind 
a substanti al amount of m olecular material. Indeed, the CO 
data of lDame et alJ (|2001 ) indicate at least three CO clouds at 
LSR velocities between and 75 kms" 1 , which presumably 
contribute to the high observed neutral column. 
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3.2. Relativistic Wind Termination Shock and Nebular 
Magnetic Field 

Knowing E for PSR J1833— 1034, we can estimate the ter- 
mination shock radius, r s , at which the pulsar wind pres- 
sure, E I {4nr 2 rjc), balances the interior pressure of the nebula, 
P„, where r] is the fractional solid angle covered by the pul- 
sar wind. Assuming nebular equipartition between particles 
and fields, P„ = B 2 J(4ir), and we have 9 S = 0."46rf ^dg 1 ^, 
where B m c is the nebular field in mG and 9, = r jd. Th e field 
derived from equipartition arguments (e.g., Pacho lcz vkl 197(1) 
for the whole PWN is B„ sa 0.3 mG (where we have used 
a s pectral index of a r = between 100 MHz and 500 GHz; 
see iBandiera et all 1200 ll and references therein), implying 
9 S w lf'67/- 1 / 2 . This can be compared with Figure [3] from 
which we infer that 9 S < 2'.'5 (note however that we do not 
actually observe a limb-brightened ring-like structure, as is 
the case in some PWNe, and it is presumably also possible 
that 9 S < 1"). The approximate match between estimated and 
inferred shock radii suggests that the nebular magnetic field 
strength is indeed relatively large, as estimated from equipar- 
tition arguments. 

For such a relatively high field, the synchrotron lifetime 
of relativistic electrons that emit 2keV photons is only ~ 
20 yr. Assuming a neb ular flow velo city profile of V(r) = 
(c/3)(r s /r) 2 (Kennel & Coroniti 1984), these electrons would 
travel ~ 12" before radiative losses become important. (The 
equivalent nebular flow ti mescale to the edge of t he P WN is 
> 500 yr.) According to Mat heson & Safi-Harbl ( 120051) . the 
X-ray spectrum steepens significantly beyond its slope close 
to the pulsar at a radius > 10"— 15" (which also appears to 
be the approximate angular scale of the relatively bright X- 
ray "plateau" seen in Fig.|2j. Additionally, it appears that the 
normalized X-ray and radio surface brightness profiles begin 
to diverge, with the nebul a becoming effectively smaller at 
X-ray energies, at > 10" ( iBock et all 120011 iDickel & Wang! 
2004, see also § 12.31 . Overall, it thus appears that a relatively 
high nebular field strength is compatible with many of the ob- 
servations. 

A break in the spectr um of G2 1.5-0.9 is inferred at 
v h > 500GHz (X h < 1 mm;lGallant & Tuffsll998t IBock et al.1 
1200 ll IBandiera et all 120011) . However, for an age > 500 yr 
and B„ w 0.25 mG, one expects a break in the spectrum due 
to synchrotron cooling at < 1 /.mi. The sub-millimeter break 
must then be due to something else, presumably reflecting the 
complex nature (or variability) of the spectrum injected by the 
pulsar or of its conversion to detectable electromagnetic radi- 
ation. 

3.3. The Evolutionary State and Age of G21.5-0.9 
For spin evolution with constant magnetic moment, the age 
of a pulsar is r = PP~ l (n - 1) _1 [1 - (P /P)"~ l ] ; the index n is 
3 for dipole braking (in th e few cases measured, 1.5 < n < 
3; e.g.. ICamilo et alJl2000l and references therein). Further 
assuming that the initial spin period Pq -C P, reduces the age 
to t c , or 4.8 kyr in the case of PSR J1833-1034. Hereafter we 
scale the age of the pulsar/SNR system by T5 = r/(5 kyr), and 
now show from multiple lines of evidence that in fact T5 <C 1 . 

The scale of structures in G2 1 .5-0.9 is R = 1 ^5^40 pc, where 
we scale angular sizes by 640 = 9/(40"). Corresponding ve- 
locities are V = 200 dsr^ 1 A840 km s" 1 , where A840 is a change 
of angular position in units of 40". For the PWN, #40 r* 1, 
and for the SNR shell, 9 40 « 150/40; thus, R PW n ~ ld 5 pc 
and Rsnr ~ 3.6afgpc. Both of these are small, similar to the 



corresponding s izes in Gl 1 .2-0.3, a 1600 year-old composite 
SNR (Kaspi e t jj .1200 1 1 ICneval ierl2005l) . 

PSR J1833— 1034 is located very near the geometrical cen- 
ter of the G21. 5-0.9 shell. Using the image presented by 
Matheson & Safi-Harb (2005), we estimate that the pulsar's 
offset from the center is A8 < 5". If the center corresponds 
to the SN site, the implied projected space velocity is Vj>sr < 
25dsT^ 1 kms" 1 . Unless t$ <C 1 , this is a very small velocity 
for a young pulsar (e.g., Lvne & Lorimen f 19941 iHobbs et alJ 
120051 according to whom the probability of Vpsr < 25 km s" 1 
is 0.45%). This reasoning was used also by Kaspi et alJ (120011) 
to argue that t <C t c for the pulsar in Gl 1.2-0.3. 

The shell in G2 1 . 5-0.9 is very faint in X-rays 
(Matheson & Safi-Harb 2005) and has not yet been de- 
tected at radio wavelengths ( Slan e et alJ 2000). This suggests 
that the system is only slightly more evolved than the Crab, 
i.e., the shell is only now becoming visible and the SNR is on 
its way to becoming a clear-cut composite like Gl 1.2-0.3 (in 
detail, such an argument requires of course an understanding 
of the circumstellar environments and explosion dynamics). 
In any case, for a uniform hydrogen ambient medium of 
density no cm" 3 , the SNR has swept up only M sw ~ 5 nod^ M Q 
and thus should be in nearly free expansion. For an ejected 
mass of M e j = \QM\qMq and SN explosion kinetic en- 
ergy of 10 51 £51 ergs, the expansion velocity would be 
3200(£ 5 i/Mio) 1/2 kms- 1 , implying r 5 = 0.25flf 5 (Mio/£ 5 i) 1/2 - 
Since a freely expanding SNR would in reality have a steep 
density profile to its ejecta, these estimates of expansion 
velocity and age are likely to be lower and upper limits, 
respectively. 

Alternatively, let us consider that the expansion has en- 
tered the adiabatic phase. Using the Sedov solution, 7?snr = 

0.31(£5i/n ) 1/5 fy/ 5 P c - If T 5 = !» thenrao/^51 = 120<i5 5 . Such 
a high ambient density is clearly inconsistent with the low lev- 
els of observed X-ray emission. If conversely we adopt typi- 
cal values of rao/Esi J$ 1, then T5 < 0.1 d 5 J 2 . But in this case 
the argument is not self-consistent, because then we would 
have M sw < M e j. The SNR is therefore still in nearly free 
expansion, for which the arguments in the above paragraph 
apply 10 . 

Consideration of PWN evolution arguments also suggests 
T5 <C 1. In general, the ratio /?pwn/^snr is initially small as 
the freely expanding ejecta move faster than the PWN expan- 
sion rate of ~ 1000 km s" 1 . However, the PWN expansion 
accelerates with the energy input by the pulsar and this ratio 
increases quickly. After the SNR ejecta decelerate and the re- 
verse shock begins to compress the PWN, the ratio decreases. 
Since the PWN in G2 1.5-0. 9 is approximately circular, is 
symmetrically placed within the SNR shell, and has the pulsar 
at its center, it is likely to be in its early pre-reverse-shock evo- 
lutionary st age. For the observed ratio /?pwn /Rrnr ~ 0.25, the 
models of Blondin. Chevalier. & Frierson (2001) are indica- 
tive of a very small age. Alternatively, considering simply 
a constant PWN expansion velocity of 1 000 V\ 000 kms -1 (for 
comparison, the Crab Nebula is expanding at ~ 1500 km s" 1 ; 
Bietenholz et al. and references therein), the observed 

PWN size implies T5 = 0.2^5^"^,. 

The PWN energetics also provide useful insight into the 

10 If M c j IOMq, as can occur, for instance, in type Ib/Ic SNe, then it 
does not follow directly from these simple arguments that M sw < M c j . How- 
eve r, in a detailed spectroscopic and spatial analysis of X-ray data on G21.5- 
0.9, Bocchino et al. 1 2005) suggest that M sw < 0.5 Mq, and that in fact the 
remnant has not yet entered the adiabatic phase. 
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age of G21. 5-0.9. According to IChevalierl j2005l) . the ratio 
Ei nt /(Er) separates "young" (ratio < 1) from older PWNe, 
where E lnt is the internal energy of the PWN. A lower limit 
on Z?jj, t can be calculated from synchrotron emission, equiva- 
lent to approximate equipartition between particles and fields, 
which is roughly con sistent with obser vations of PWNe. We 
use equation (46) of Chevalier (2005) to calculate that for 

G2 1.5-0.9, E mia ~ 2 9 x 10 47 t/f /7 erg s. Here we have used 
Pi = 1 +2a r » 1 04 dSalter et al.lH989h and p 2 = 1 +2a x « 3 
(ISafi-HarbetalJl2001l) . as well as u b = 500 GHz and L vb w 
1.2 x 1 23 c/f ergs s" 1 Hz" 1 , where p\ and p 2 are the energy in- 
dices for radio- and X-ray-generating particles, respectively, 
and L v b is the spectral luminosity at v b . With the measured 
E for PSR J1833-1034, we have Et = 5.3 x 10 48 r 5 ergs, 
which is larger than Ei nt (typically within a factor of a few 
ofEmin) for any reasonable age. This then implies r <C t c (see 
IChevalierl2005T) . 

These various arguments all lead to the conclusion that 
for PSR J1833-1034 and G21.5-0.9, r 5 < 1, and maybe 
r 5 lOOOyr (a similar co nclusion was reached independently 
bv IBocchino et al.ll2005l who in an analysis of X-ray data 
on G21.5-0.9 obtain an age range of 200-1000yr). In 
turn this suggests that Po > 55 ms essentially independent of 
braking index. Measuring n would of course be very in- 
teresting and would further constrain the rotational history 
of J1833-1034, but this is often impossible. For example, 
for PSR J0205+6449 in SNR 3C 58, the rotational insta- 
bility of the neutron star prevents n from being determined 
( Ranso m et aTl2"0 04 ) . 

While the age inferred here would make PSR J1833-1034 
among the 2-3 youngest known neutron stars in the Galaxy, 
it is highly unlikely that its progenitor SN was observed as 
an historical event (and ind eed no SN is recorded in t his di- 
rection in the compilation of Stephenson & Green 2002). The 
observed hydrogen column density corresponds to a visual ex- 
tinction of 10-11 magnitudes. The peak visual magnitude of 
the supernova might thus have been only 5-6 (for an absolute 
magnitude My = -18.5), which is likely to have gone unno- 
ticed. 

3.4. Comparison of G21. 5-0.9, 3C 58 and their Pulsars 

The measured spin parameters of P SR Jl 833-1034 are very 
similar to those of PSR J0205+6449 (Murra v~l.l2002l) . the 
central engine in 3C 58 (another famous SNR long-classified 
as Crab-like, from which faint thermal emission has no w been 
detected; see IBocchino etalJl2001t ISlane et all 12004). It is 
therefore of interest to compare the properties of the respec- 
tive PWNe, and we list some important parameters of both 
systems in Table [3] 

Most noticeably, 3C 58 is a much larger and more asymmet- 
ric PWN than G2 1.5-0. 9, and its ratio of X-ray luminosity to 
spin-down luminosity, Lx/E, is ss 11 times smaller. The in- 
ferred pulsar wind shock radii differ by a factor of > 3, but this 
is understandable if the nebular pressures differ by an order of 
magnitude, which may well be the case, as inferr ed fr om the 
respective equipartition field estimates (see also § I3.2> . Also, 
the break frequency v b for 3C 58 seems to be ~ 10 times lower 
than for G21.5-0.9 11 . 

" This may or may not be of significance. Spectral indices have not been 
reliably measured just below and above the putative v\, for either of these 
systems, and a "break" frequency is instead inferred from extrapolations of 
fluxes measured 1-2 orders of magnitude apart. The complexity of the broad 
band emission from these PWNe may not be encapsulated by references to, 



TABLE 3 

Parameters of PWNe G21.5-0.9, 3C 58, and their Pulsars. 



Parameter G21.5-0.9 3C58 



Distance, d (kpc) 4.7 3.2 

PWN size (pc 2 ) a 1.8x1.8 10x5 

Termination shock radius, r s (pc) < 0.05 m 0.15 

Radio luminosity (10 7 -10 u Hz), (10 34 ergs s -1 ) . . L4 3.0 

Radio spectral index, a r 0.0 0.1 

Break frequency, v h (GHz) > 500? ~ 50 

Minimum internal energy, £* m ; n (ergs) ~ 2.5 X 10 47 ~ 10 48 

X-ray luminosity (0.5-10 keV), L x (10 35 ergs s" 1 ) 2.8 0.2 

X-ray spectral index, cex b 0.4-1.3 0.6-1.6 

Equipartition magnetic field strength, B„ (mG) . . 0.3 ~ 0.08 

Age, r (kyr) < 1000? 824? 

Pulsar period, P (ms) 61.8 65.6 

Spin-down luminosity, E (10 37 ergs s~') 3.3 2.7 

Characteristic age, r c (kyr) 4.8 5.4 

Initial spin period, Po (ms) > 55? ~ 60? 



NOTE. — We use the nominal distances listed here to infer all quantities 
that depend thereon. Parameters for G2 1.5-0. 9 are from this work and its refer- 
ences, while those for 3C 58 are from Slane et al. 1 2004) and references therein. 

a At least in projection, G21.5-0.9 is approximately round (Fig.[2j and 3C 58 
is elongated. In both PWNe, the X-ray nebula is so mewhat smaller tha n the 
radio nebula, an effect more pronounced in G21.5-0.9 i Bock et al. 2001 ). 

b Both X-ray spectra steepen radi ally from the center i n sim ilar fashion, 
within the ranges of a x ind icated <Safi-Harb et all tZOOl: Slane et al. 2004; 
Matheson & Safi-Harb 2QQ3). 



3C 58 is widely thought to be the remnant of SN 1 181 CE 
( Steph enson & Greenl 120021) . In order to bring the his- 
torical age and m e asured spin parameters into agree- 
ment, iMurrav et alJ J2002I) proposed that the initial pe- 
riod of the pulsar was Po ~ 60 ms , co mparable to what 
we infer for PSR J1833-1034 (§ E3- On the other 
hand, good arguments suggest that 3C 58 may be older 
jBietenholz. Kassim. & Weilerl 1200 ll) . in which case Pq <C 
60 ms. In any case, different spin histories cannot explain 
different present-day values of Lx, since the relevant elec- 
trons/positrons cool relatively quickly in both PWNe. 

If its apparently lower inferred v b and steeper X-ray spec- 
trum are of significance, much of the pulsar power in 3C 58 
may be deposited at lower particle energies, and hence radi- 
ated at lower frequencies, compared to the case for G2 1 .5-0.9. 
In addition, the lower nebular magnetic field strength inferred 
for 3C 58 (see Table[3} means that its relativistic particles are 
relatively poor radiators, and this may explain the vastly dif- 
ferent X-ray efficiencies in the two systems. 

But why is the nebular field strength in 3C 58 apparently 
so much lower than for G2 1.5-0. 9? If 3C 58 is not the rem- 
nant of SN 1181 CE bu t is substantially old er, as argued by 
iBietenholz et alJ J2001I) and lChevalieil <2005l) . then the larger 
PWN size would follow naturally, and its field would also 
have decayed. Alternatively, even if the ages for both systems 
are similar, the PWN radius at early stages depends on SN 
and pulsar properties as R oc £ 1/5 £'sn 3 W M e f 1 ' 2 , where Esn 
and M e j are the kinetic energy and eje cted mass, respectively, 
of the associated supernova explos ion iRevnolds & ChevalieJ 
Il984tfvan der Swaluw et all20 01). Therefore, it may be pos- 
sible to obtain a larger size and lower magnetic field for 
3C 58 if either the supernova explosion was more energetic 
or the ejected mass smaller than the case for G2 1.5-0. 9, ei- 
ther of which are reasonable. Additionally, perhaps insta- 

e.g., a single value of v/,. 
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bilities in the magnetic field downstream of the termination 
shock are also of relevance for explaining the elongated shape 
of 3C 58 and maybe its size (see, e.g., van der Swaluw 2003; 
Slane et a l 120041) . 

Finally, we note that the order of magnitude discrepancy 
in L x /E between PSRs J1833-1034 and J0205+6449, two 
pulsars with near-identical spin parameters, makes it clear 
that the predictive power of simpl e relations between E 
and L x for pulsars and their PWNe JSeward & Wang||1988l 
Bec ker & Trumperll997t iPossenti et alJ2002l) is limited. 

3.5. The Radio Spectrum of PSR J1833-1034 and Young 
Pulsars 

The radio spectrum of PSR J1833-1034 between 0.8 GHz 
and 2.0GHz is un usually steep (a « 3.0; see § 12. 2\ . 
lLorimer et al.l (Q995) show that the mean spectral index for 
pulsars is w 1.6, with a suggestion that it is smaller (~ 1) 
for young pulsars. The Crab pulsar is an exception, with 
a = 3.1. PSR J1833-1034 has a pulse duty cycle of 4%, 
smaller than virtually all pulsars with P < 0.1 s (but simi- 
lar to that of J0205+644 9 in 3C 58, which h as similar pe- 
riod as well, and a ~ 2; Camilo et al. 2002d). With such a 
steep spectrum, the pulsar is easier to detect at low frequen- 
cies, such as 820 MHz at GBT However, for good reasons, 
the first frequency of choice to employ in searching for dis- 
tant pulsars along the Galactic plane remains 1400 MHz. At 
this frequency, PSR J1833— 1034 has a very low flux density 
and its luminosity L1400 = Si4oo<^ 2 ~ 1 -8 d\ mJy kpc 2 is among 
the smallest known. The detection of such pulsars is therefore 
not easy. Nevertheless, this may owe more to the state of our 
detection technology than to the intrinsic rarity of such "faint" 
neutron stars. 

We have compiled information about all known young 
rotation-powered pulsars, defined here somewha t arbitrarily 
as tho se with r c < 30 kyr (see the catalog of Manchest er et al.l 
120051 but some of the distances used here are from orig- 
inal references). There are 36 such pulsars, of which 33 
have been detected at radio wavelengths. Of these, 23 
have a "high luminosity" L1400 ( 1 8-140 mJy kpc 2 ; median 
56mJykpc 2 ): 14 have been long known (of which three 
were first detected as X-ray pulsars) and nine were discov- 
ered in the Parkes multibeam Galactic plane survey (e.g., 
iKramer et al.l l2003b). With PSR J1833-1034, we know to- 
day of 10 young "low-luminosity" pulsars (0.5-6.3 mJy kpc 2 ; 
median 2.4mJykpc 2 ). Seven of these were detected in deep 
searches of PWNe Edpern et al.ll2001j: iMurrav et alJl2002L 
ICamilo et al J2002alblcldl Roberts et al. 2002} and this work), 
two more were detected in the Parkes multibeam survey 
(Camilo et al. 2001, 2004), and only one has been long known 
(B1853+01, discovered in a direct ed search of SNR W44; 
Wolsz czan. Cordes. & Dewe\ll991l) . Evidently, such pulsars 
are not so rare, but merely require more effort to uncover. 

4. SUMMARY 

We have discovered the 61.8 ms pulsar J183 3— 1034 at the 
center of the composite SNR G21.5-0.9 (§12. 1>. Second only 
to the Crab in spin-down luminosity among known rotation- 
powered Galactic pulsars, with £ = 3.3 x 10 37 ergs s" 1 (§ 12. 2> . 
PSR J1833— 1034 was born in the SN explosion that gave rise 
to G2 1.5-0.9, which we estimate from the properties of the 
pulsar, PWN, and SNR shell occurred < 1000 years ago. This 
is substantially smaller than the pulsar characteristic age and 
in the simplest interpretation implies a birth spin period of > 



55 ms (§ 13. 3> . providing further evidence for a relatively large 
range in initial periods for neutron stars (e.g., Migliazzo et al. 
20021 IKramer et all2003al) . 

PSR J1833— 1034 is exceedingly weak at the radio fre- 
quency of 1.4 GHz used at the Parkes telescope to first de- 
tect it. We have revisited the distance to G2 1.5-0. 9 based 
on H I and CO obse rvations, and find the best estimate to be 
4.7 ± 0.4 kpc (§ 13. 1> . for which the pulse-averaged luminosity 
at 1 .4 GHz is about 1 .5 mJy kpc 2 , one of the smallest known 
for a young puls ar. Observations at other frequencies with 
the GBT (§E3 suggest that J1833-1034 has an unusually 
steep spectrum, both of which findings have significant im- 
plications for the detection of additional very young pulsars 
and for studies of the radio lumin osity , beaming fraction, and 
birth rate of such neutron stars (§ 13. 5> . 

We have analyzed a large number of Chandra X-ray obser- 
vations of the system and obtained the best image to date of 
the inner core of G2 1.5-0. 9. This shows clear evidence for 
a pulsar component offset from the center of a surrounding 
very compact elliptical emission region t hat l ikely bounds the 
termination shock of the pulsar wind (§ I2.3> . From this, un- 
der the assumption of energy equipartition in the PWN, we 
have obtained an estimate of the nebular magnetic field that 
is consistent with other determinations and plausibly is also 
broadly compatible with the X-ray and radio surface bright- 
ne ss an d spectral energy distribution observed in the PWN 
(§ 13.21 . PSR J1833— 1034 has present-day spin parameters 
very similar to those of PSR J0205+6449 in SNR 3C 58, 
while the respec tive PWNe exhibit some significantly differ- 
ent properties (§ I3.4> . 

We have searched fo r X-r ay pulsations from PSR J 1833- 
1034 without success (§ 12.41 . Our pulsed fraction limit for an 
assumed sinusoidal pulse shape is about 70% of the point-like 
luminosity of 3 x 1Q~ 5 E. Depending on the rotational stabil- 
ity of the pulsar, to be determined from timing observations, 
we may in future obtain a better limit or, if the X-ray beam 
intersects our line of sight, a detection of such emission. In- 
terestingly, no EGR ET 7-rav source is known at the location 
of J1833-1034 i lHartman et alii 19991) — despite the fact that 
this neutron star has the fourth largest spin-down flux (E / d 2 ) 
among known rotation-powered pulsars (after the Crab, Vela, 
and J0 205+6449, which also is not a known 7-ray emitter; 
Cami lo et alJl2002dl) . Assuming it beams toward the Earth, 
its efficiency for converting rotational power into > 100 Me V 
7 rays is relatively low, in keeping with a suggested inverse re- 
lation with E llThompson et alll999l) . Certainly, PSR J1833- 
1 034 should be a prime target for the future Gamma-ray Large 
Area Space Telescope. 
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Phase Phase Phase 

FIG. 1. — Radio pulse profiles of PSR JI 833— 1034. Central frequencies increase from left to right: 820 MHz (GBT), 1374MHz (Parkes; discovery observation), 
and 1850 MHz (GBT). In each plot, the linear grey-scale represents pulsar strength as a function of rotational phase and time (increasing from bottom), with the 
total observation split into and displayed as 64 sub-integrations. At top is shown the summed pulse profile for the full observation. Two pulse periods are shown 
in all cases. The 1850 MHz profile was obtained from a bandwidth of 400 MHz (the pulsar was not detected in the upper 200 MHz of bandwidth recorded with 
SPIGOT; see TablelTT. and the white band indicates removal of data due to interference. There is evidence of apparent interstellar scintillation in the time- varying 
flux densities of the 820 MHz and possibly 1374 MHz data sets (and the nominal flux densities for separate 820 MHz and 1374 MHz observations differ by 
~ 10%— 30%; see TablelTT. The observed pulse full-width at half-maximum is 3.0ms (0.05P) at 1374MHz and 2.5 ms at the other two frequencies. Smearing 
due to differential dispersive propagation within one frequency channel is 1.3 ms for the 820 MHz data and 1.6 ms for the 1374 MHz data (which, additionally, 
have a sample interval of 0.8 ms). Therefore, the intrinsic pulse FWHM is very nearly 0.04P at all frequencies. 
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FIG. 2. — False-color Chandra ACIS-S image of the PWN in SNR G2 1.5-0. 9 based on 58 ks of data, showing the elongated core embedded in the more diffuse 
nebula. The SNR shell lies beyond the edge of this image. The image has been smoothed by two pixels (0"98), and the intensity levels increase logarithmically 
from a surface brightness of 9ctsarcsec~ 2 in the outer regions (blue) to a value of 1640cts arcsec~ 2 in the core (white). See ij |2.3l for further details about the 
image. 
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FIG. 3. — Grey-scale Chandra HRC image of the central core in G21.5-0.9 based on 351 ks of data from multiple observations (compare to Fig.[2j. The image 
has been smoothed by two pixels (0"26), and any residual registration error/blurring is at the level of < 0"3. A compact source, presumably corresponding to 
PSR J 1833-1034, is surrounded by an elongated emission region. The grey-scale intensity increases linearly from 520 cts arcsec~ 2 at the outer edge of the core 
(e.g., near R.A. 18 h 33 m 33?5, Decl. -10°34'05"5) to 2500 cts arcsec~ 2 at the position of the compact source. See §j}|23landl2~4lfor more details. 



